DOCOHENT BESOME 



BD 195 U32 SE 033 593 

AOTHOR Posemaiir Leonard D. 

TITLE Light and Sound: Evolutionary Aspects. Physical 

Processes in Terrestrial and Aquatic Ecosystems^ 
Transport Processes. 

TNSTITOTION Washington Univ., Seattle. Center for Quantitative 

Science in Forestry, Fisheries and Wildlife. 

SPONS AGENCY National Science Foundation, Washington, D. C. 

POB DATE Jan 78 

GPANT NSF-GZ-2930: NSF-SED74-1 7696 

NOTE 93p.: For related documents, see SE 033 581-597. Not 

available in hard copy due to marginal legibility of 
original document. 



MF01 Plus Postage. PC Not Available from EDRS. 
Acoustics: *Bioloqy: college Science: Computer 
Assisted Instruction: Computer Programs : Ecology: 
Energy: Environmental Education: ^Environmental 
Influences: *Evolution: Higher Education; 
Instructional Materials: *Interdisciplinary Approach: 
Light : Physical Sciences: *?adiation ; Science 
Education: Science Instruction 

ABSTRACT 

These materials were desianed to be used by life 
science students for instruction in the application of physical 
theory tc ecosystem operation. Host modules contain computer programs 
which are built around a particular application of a physical 
process^ This module is concerned with the exchange of energy between 
an organism and its environment in the form of radiation. Classical 
and modern radiation theory are discussed and the applications of 
these physical principles to analysis of biological systems are 
presented. E'lrphasis in these applications is upon the evolutionary 
significance of the physical processes, in one type of application, 
physical theory is used to isolate constant factors in the 
environment tc which all orqanisms must adjust their evolution, in a 
second application, physical theory is used to elucidate the 
constraints governing the evolution of biological systems, it is 
demonstrated that radiation theory is applicable to a wide range of 
sense organs, radiation types, wavelengths, and organisms. 
(Author/CS) 
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INTRODUCTION 



All living organisms interact with their environment through a 
continuous exchange of matter and energy. Conceptually, this can be 
thought of as a consequence of the "openness" of the system in a thermo- 
dynamic sense. For a complete understanding of the biology of an 
organism, it is essential to determine the ways in which the organism 
interacts with the biotic and abiotic factors present in the external 
enviroiiment . To most biologists > the exchange of energy between an 
organi/?.;.! and its environment in the form of radiation is probably the 
least well understood abiotic factor affecting the organism. As D. Gates 
has pointed out,^' this exchange is also the most difficult to measure. 
However, its importance to the organism cannot be overemphasized. 

There are two ways in v/hich the process of radiation exchange is 
of biological interest. The governing principles of this physical 
process are different for the two distinct biological applications. 

First, there is the dynamical balance of the radiation flow between 
the organism and its environment, and its effect on the homeostasis of 
the organism. This dynamical balance is best described by the modern 
form of radiation theory, which is based on Planck's theory of black 
body radiation, and which is exemplified in the generally useful results 



[Gates, D. M. "Radiation incident on an organism" 1978, p. 1] 
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entailed in Wieii^s law of shift and the Stef an-BolfcEinann law. Several 
simple biological applications of this modern radiation theory are 
developed in this module. Further applications are present . j •^^j'^where 
in this series (Gates 1977,1978; Hatheway 1977, Stevenson; 1/' ' . .. 

Second, there is the interpretation of the energy he 
organism, as a means of obtaining information about the Gnv^r.. 
The selective advantage that results from improved information i_ - ...aring 
has led to the evolution of many specialized sense organs. The function 
of these organs is best described by the classical theory of ra'^lation, 
particularly as entailed in Rayleigh^s criterion and the Doppler affect. 
The balance of biological examples discussed here are of this second 
type. 

The nature of the physical processes of radiation exchange imposes 
strong constraints on the evolution of biological systems. This module 
presents some of the results of classical and modern radiation theory and 
shows how they elucidate these physical constraints. It is demonstrated 
that radiation theory is applicable to a \vide range of sense organs 
(eyes, ears, pits); radiation types (electromagnetic, sound), wavelengths 
(optical, infrared), and organisms. 

BLACK BODY RADIATION 

In any discussion of radiative energy exchange, it is necessary to begin 
v;ith an idealization of this process, ihich is described by the theory of 
black body radiation. This is rather awkward in the present instance 
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since the biological. implications of this abstraction from the real 
world are not inunediately obvious. In fact, some of the most important 
aspects of biological systems result from their deviations from the 
idealizations of physics. For example, most organisms are not perfectly 
''black" in the physical sense of the word; that is, they do not absorb 
and emit radiation perfectly at all frequencies. The effect of the 
coloration of organisms will be to vary the degree to which the black body 
idealization can describe the systerr. However, it is the nature of 
science to explain in detail the simple system in the hope that the more 
- complicated systems can be explained by small variations of the simple 
theory, and the expectation that the major results will be carried over 
to the more complicated systems. For radiative energy exchange this is 
precisely the case, and thus there is sufficient justification for present- 
ing a discussion of black body radiation which, though not "simple" in the 
comon sense of the word, is essential for an understanding of the more 
complicated systems which are the substance of biology. 

The theory of black body radiation had a lengthy and fascinating 
development during the ISOO*^-, culminating in the work of Planck at the 
turn of the century. Many of the results discussed below were known an 
various stages during this period of repeated empirical research and 
hypothesis testing. However, the important results for biological 
applications can be derived directly from the formula of Planck, which 
mathematically (though not physically) culminated this field of experi-- 
mental physics. The physical explanation of Planck's formula awaited 

e 



the work of Einstein, and ultimately the early founders of Quantum 
Mechanics . 

Planck's problem was to specify the energy distribution of the 
elfictromagnetic radiation within a cavity, the walls of which are at 
thermal equilibrium at the temperature T. That is, the problem is to 
find the energy density E(v) (J s m"^) of " the radiation within the 
cavity between the frequencies v and v + dv. The shape of the curve E(v) 
vs. V had been determined experimentally at various temperatures 
[Fig. 1 ] . Planck found a solution to this problem in the form: 

E(v) = Slllji p/kT _ ^ 

where v is the frequency of the radiation, T is the absolute temperature 
of the walls of the cavity, and h, k and c are physical constants [see 
Appendix IJ. This energy distribution is called the black body energy 
distribution because it is the same as the energy distribution of radiation 
emitted by a perfect.'y black object which is at absolute temperature T. 
For measurement purposes, it Is more useful to express this energy density 
in terms of the wavelength A. 
Problem 1 

Show that the energy density per unit wavelength is given by: 

STThc , hc/XkT "1 
E(A) = — [e - 1] 

. . Hint: A = c/v 
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FIGURE 1» BLflCKBODY RRDIRTION PLRNCKtS LRW VS. OBSERVED 
DRTR FROM LUMMER RND PRINDSHEIIi (1901) 

2500r 




WRVELENGTH IN MICRONS 
CURVES- PLRNCKtS LRW FOR T:= 1646, 1449, 1259, 1095, 998, 904, 723 
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Welti's Law of Shift 

It can be shown that the E(X) distribution can be described by the 
wavelength at which the distribution reaches its maximum. Qualitatively 
increasing the temperature T shifts the maximuia of the E(X)[E(v)] curve 
to a smaller X (higher v). 
Example 1 

The wavelength at vjhich the energy density is a maximum at a given 
temperature T can be found as follows: 
From problem 1: 

-1 



E(X) = 



Srrhc 



hcAkT 
e - 1 



Make the substitution u = hc/XkT 

5 ^ 



Ea) = 



SjrkV 
c h 



e" - 1 



From calculus, the maximum can be found by solving: 



4^ = 0, which gives 
du 



dE STTk^T^ 



du 



c h 



(5u ) (e -1) - e (u ) (e -1) 



= 0 



or 



(5u^e"-l) - uV)(e"-l) 

j.Au 4 5u ^ 
5ue -5u - ue = 0 



-2 



0 which implies 
or finally: 



-u 



e + l/5u -1=0. 



This transcendental equation can bG solved by numerical techniques (see 
Appendix 2) to give: 

u = 4.9651 
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Thus 

X T = = 2.8978 X io"\ K « b. 
m uk 

b is called the Wien constant and the equation 

X T = b 

m 

is called Wien's law of shift (1896 Wilhelm Wien). Therefore, the 
maximum of the black body curve E(X) at temperatures T^^T^yTy... falls 
at . . . so that 

Vl = V2 = V3= ' 
[see Fig. 2]. 

The Wein lav/ of shift is a very powerful result because it embodies 
a single relation betv/een wavelength and absolute temperature. This is 
exemplified in the folloicing problem. 
Problem 2 



Assuming the following objects radiate as black bodies, what is the 



\^ for: 



a) the sun at T = 5700 K 

s 

b) a mammal with a surface temperature T - 37°C 

c) a lizard with a surface temperature T = lO^^C 

d) the clear night sky, radiation temperature T - '-33°C 



Answers; 



a) 5.08 X lO^'^m = 5080 A (yellow) 

b) 9.66 X 10"\i (infrared) 

c) 1.024 X 10~\ (infrared) 

d) 1.21 X 10"\ (infrared) |0 

'} r 
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Stefan-Boltzmann Law 



-3. 



If E is the total energy density in all frequencies (J m ): 

Sirh 



E = E(v)dv 
Making a change of variable. 



3 

V dv 



hv/kT 
o e - 1 



kT 



CD 



X = hv /kT which implies dv = — dx. 

h . 



Equation (1) becomes 



E = 



kT 



3^ 
X dx 

e -1 



(2), 



The value of the integral in Equation (2) 



I = 



3^ 
X dx 

o e -1 



(3) 



can be obtained by numerical integration. This is a technique for 
approximating a definite integral by a finite sum. The value of the 
integral in Equation (3) ir> found in Appendix 2 to be: 

I = 6.4983. 
Equation (2) can be I'e-expressed as 

4 



E = 



87Th ,k. 
3 ^h^ ^. 



4 4 

T = aT 



(4) 



where the constant a is equal to the value of tlie coni^Lants in the 
square brackets in Equation (4). The value of a is found to be 

a = 4 (6.4983) = 7.5-643 x 10"^° J m'^ k'^ 

c h 

The energy radiated by a black body per unit area per unit time. 



that is, the radiation emittance E , can be sho^^m to be 

n 
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E = OT -(5) 
n 

where E = radiation emittance and 
n 

a = (^)ca = 5.6693 x lO^^W m^^K"^. 
Equation (5) is called the Stef an-Boltzmann law and the constant O is 
called the Stef an-Boltzraann constant. Although Equation (5) is quantit- 
atively precise, the functional relationvship of the variables is more 
important than the exact numerical value of the constants. Of particular 
interest is that the radiation emittance is proportional to the fourth 
power of the absolute temperature. 

The equation is of direct biological interest when it is realized 
that the Stef an-Boltzmann law governs the radiative energy exchange 
between an organism and its environment. Although this ^exchange of energy 
is usually not considered when calculating the energy budget of an organism 
it can in some circumstances be a major factor affecting the heat balance 
of an organism. The following example and problem constitute a brief 
introduction %o the application of this physical law to biology. For a 
more detailed exploration of this relationship, see Gates (1978) and 
Hatheway (1977). 
Examj)le 2 

A bat with a surface temperature of 37**C is living In a cave which 

has walls at a temperature of 12^ C. I^at is the radiation emittance of 

2 

the bat? If the bat has a total surface area of 100 cm , what i« its rate 
of heat loss by this process of radiative exchange? 
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Solution: Define 

S the bat's surface temperature (in K) 

T = the ambient temperature (i.e., the temperature of the cave) 
a 

(in K). 

Clearly the bat is emitting radiant energy at the rate: 
4 

E = aT 
o o 

and receiving radiant energy (from the walls) at the rate: 

E = qT ^. 

a a 

So the overall energy exchange by this process is given by: 

E = a(T - T ) 
o a 

with the convention that E is positive when the overall radiation balance 

is such that the bat is emitting more than it absorbs; i.e., it is losing 

energy to the environment at the rate E. 

For this example: 

T = (273 + 37) K = 310 K 
o 

T = (273 + 12) K = 285 K 
a 

which gives an energy emittance of: 

E = 5.67 X 10'"^(W Di"^ K"^) [(300 K)^ - (285 K)^] 
^ 5.67 X 10"^[9.235 x 10^ - 6.598 x 10^] W m ^ 
= 5.67 [92.35 - 65.93] W rrT^ 
= 149. 5w m'^ 

The bat's total heat loss rate H, by this process is 
H = E(SA) 

where SA = surface area, so 

-2 2 
H 149.5 W m (100 cm ) 



= 1.495 W. 



12 



Although in this example the heat loss rate was relatively small, 
when there are larger differences between the organism's temperature and 
the (radiation) temperature of the environment, the heat loss rate can 
be quite high. This situation pertains in the following problem. 
Problem 3 

a) The horns of a reindeer are at approximately internal body 

temperature, though the rest of the deer's body is well insulated. If 

these horns are at a temperature (T^) of 32°C, and they have a surface 
2 

area of 1000 cm , what is the heat loss rate if the environment is at a 
(radiation) temperature (T^) of -10°C? Convert your answer to calories/ 
minute. 

[See Gates 1968 for some actual organismal temperature measurements.] 
Answer : H:^ 22 W, 

b) A bald mountaineer is outside on a clear night. Although he is 

heavily bundled in his down parka, he has forgotten his hat. If his bald 

2 

pate has a surface area of 75 cm and is at a temperature of 32''C, what 
is his heat loss rate to the night sky which is at a (radiation) temper-- 
ature of -45°C? [Is this why they say, "if your feet are cold, pat on 
your hat."?] Convert to calories/minute. 
Answer: Yl-^ 25.3 W = 362 cal/min. 
The solar constant 

A physical parameter which imposes a significant constraint upon 
the evolution of terrestrial organisms is the amount of energy generated 
by the sun which is intercepted at the earth's surface. Almost all of 

i 

the energy which drives the biological processes on the earth comes 



is 



13 



directly or indirectly from the sun. Thus the amount of the sun's energy 
that the earth receives is of interest for at least two reasons: (1) This 
is the upper limit to the amount of energy that could be absorbed by a 
surface and turned into heat. (2) This is the upper limit to the amount 
of energy that could be converted by a photosynthetic unit into chemical 
energy. Since the actual amount of sunlight available at the earth's 
surface is dependent upon absorption by the atmosphere, and thus on such 
variables as latitude, the sun's altitude in the sky (see Gates 1978, 
section on Direct Solar Radiation), or cloudiness, the measure of the 
sun's energy available at the earth is determined instead for a unit area 
just outside the earth's atmosphere. This measure of the sun's energy 
(actually the rate of the available energy, i.e., power) is called the 
solar constant . 

The solar constant can be calculated theoretically from the 

Stefan-Boltzmann law and conservation of energy. Assuming the sun to be 

a black body with a surface temperature T^, by the Stef an-Boltzmann law, 

its radiation emittance (energy radiated per unit area per unit time) is 
given by: 

^s^^^s''- (5) 
Thus the total power emitted by the sun can be found by integrating E 

over the sun's surface. Since T^ is (assumed to be) constant, this is 

equivalent to multiplying the constant by the surface area of a sphere 

with the sun's rad ius (r ). The total power is given by: 

E = 47irJ E^. ^ (6) 

By conservation of energy, if no energy is lost in the intervening space, 

all of the radiation emitted r^adially out from the sun will be intercepted 

■ 
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by an (imaginary) sphere centered on the sun with a radius equal to the 
mean distance from the sun to the earth (r_) . The power flow per unit 
area (E_) across this (imaginary) sphere will be equal to the total 

2 

power available (E) . divided by the surface area of the sphere (Airr-.). 



That is. 



Eg = E/4Trr^. * (7) 



Substituting Equation (6) in Equation (7) gives 
4:rr E r 

The quantity E which is the power available (per unit area, per unit 

time) on a surface at the earth *s distance from the sun, is the parameter^ 

of interest, the solar constant. 

Inserting Equation (5) into Equation (8) gives immediately: 
2 



r 

s 

2 



aT ^. (9) 

S 



Problem 4 

Use Equation (9) and the following physical parameters to calculate 

the 5;olar constant: 

T = 5700 K 
s 

r = 7.1 X lO^m 
s 

r-, = 1.49 X IQ-^-^m 

a = 5.6693 x lo""^ W m"*^ k""^ 

-2 . -1 

Convert your answer to cal cm mm 

3 -2 -1 -2 -1 
Answer: E_ = 1.36 x 10 Win. s <Dr 1.95 cal cm min *. 
— — t, 
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III. Resolution 

For most problems of biological interest (e.g. the shape of the 

lens of the eye), only the particulate nature of light, especially its 

rectilinear propagation, is important. This can be thought of as a 

result of the large size of the elements of biological interest compared 

to the wavelength of light. For example, for the sun: A 5 x lo""^m 

m 

[problem 2a], which is much smaller than the physical dimension of any 
organism, or even any organ in any organism. Only in certain specialized 
biological research, for example the investigation of cell ultrastructure, 
does the wavelength of light become a limiting factor. As the dimensions 
of the cellular constituents approximate that of light, instruments 
capable of observing at shorter wavelengths (e.g. electron microscopes) 
are required. 

There is at least one other area of biology where the physical 
dimensions of the objects of interest are small enough that the wave 
nature of light becomes an important consideration: the resolution of 
incoming radiation by sense organs. 

The term resolution is used here to denote the ability to 
distinguish two objects. That is, how close together must two objects 
be before they appear as a single object to the organism? Clearly, 
this depends to some extent on the organism in question; eagles resolve 
better than do moles. Also, the limit of resolution vjill be affected 
by the experience of the observer. Nevertheless, it is useful to have 
some rough comparison of the resolving ability of various systems, and 
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for this purpose Rayleigh's criterion is most often employed. 

Rayleigh's criterion is derived from the theory of Fresnel diffract- 
ion by a circular opening (for a derivation- and detailed discussion, see 
Elmore & Heald 1969). Qualitatively, this criterion is seen to be a 
measure of the relative overlap of the diffraction patterns of the images 
of objects on the focal plane of the observer. More specifically, it is 
a measure of the displacement of the central maxima of these diffraction 
•patterns. Two images are said to be just resolved when the central 
maximum of one diffraction pattern falls on the first minimum of the 
other (see Fig. 3) . 

From Figure 3, it can be seen that the limit of resolution is some 
distance fix between the two images on the focal plane, or alternatively 
some angle 66 between the central axis of the lens and the ray connect- 
ing the center point of the lens to the centers of the principal maxima 
on the focal plane [pj[,P2 of Fig. 3]. Therefore, in comparing the 
resolving power of different instruments the reciprocal of these quantities 
should be used (i.e., l/6x or 1/66), so that the 'faore powerful" instrument 
is able to distinguish objects which are unresolved by the "less powerful" 
instrument. For a distant object, the objective (lens) is the aperture 
which creates the circular diffraction. The mathematical details of 
Fresnel diffraction lead to the result that the position of the first 
minimum relative to the principal maximum will be: 

sin a = 0.61 X/r (11) 

where 
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Figure 3: The Rayleigh criterion: Fresnel diffraction by circular 
aperture. 




a) A simple diffraction apparatus. The lens L fornfe^the images of the point i 

sources and on the focal plane F. These images are not points but ^ 
disks centered at Pj^ and P^ whose intensity is greatest at the center as 

indicated by the "humps" behind the focal plane. [From Sears and Zemansky p. 619 



b) 




Diffraction patterns almost coincident: unresolved . 



6x 



/ 




c) Diffraction patterns v/ith the first niiniiauni of the image of object 2 falling on 
the principal maximum of the image of object 1: just resolved . (This is the 
Rayleigh limit.) 
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X - the wavelength of radiation which forms the image, 

r = the radius of the objective, 

a ^ the angle between the images of the objects. 

The limit to resolution given by Eq. (12) is known as Rayleigh^s 
criterion. Using the familiar (Taylor) expansion of the sine function 
one finds that: 

3 5 7 
sm a = a--a +a -a +••• 

3 

or sin a = a + 0(a ), 

Here 0(a ) means terms of order (size) strictly less than a^. If 
3 

a IS small, then a is very small and one can make the small angle 
approximation 

sin a <c: a 
which when substituted into Eq, (11) yields: 

a = 0.61 X/r , (12) 
Rayleigh^s criterion in the simple form of Eq, (12) is a powerful 
tool in an analysis by the biologist of the evolution of sense organs. 
In particular, such questions as the scaling of eye size with organismal 
s^ze (e.g., why are a whalers eyes the same size as those of a cow?), 
or the cellular organization of the organs themselves (e.g* the distri- 
bution of cones on the retina) have relatively simple answers when 
analyzed in this way. The ramifications of this analysis are explored 
in the following examples and problems. 
Example 4 - 

From the Rayleigh criterion, the results of problem 2a, and the 
fact that in bright daylight the diameter of the human pupil is*<t'2.0 mm. 
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determine the resolution of the human eye in minutes of arc. What does 
this resolution limit correspond to at 10 m? 
Solution: A^^ = 5.08 x lO^'m (from prob. 2a) 
r = 1.0 mm = 1.0 X 10 m. 
Putting these parameters into the formula for the Rayleigh criterion, 
Eq. (12), v/e get: 

a = .61 A/r 

= 0.61(5.08 X 10"^m)/(1.0 x lo""\) 
-3 

- 3.1 X 10 radians. 

Since 27t radians = 360° = 2.16 x 10^ min, we know that 

3 

1 radian <^ 3.43 x 10 min. 
So a = (3.1 X 10^ min/rad)(3.1 x lo""*^ rad) 

min of arc. 

At 10m this angle will correspond to an object separation of: 

a = S /d 
o 

where S = separation 

o ^ 

d ~ distance 
or = ad = (3.17 X lo"^) (10m) = 3.117 x io"\ 
- 3.17 mm. 

Thus, in daylight the human eye can- just barely resolve tv/o objects which 
are 3.17 mm apart from one another and 10m distant. 

It should be noted at this point that the Rayleigh criterion is a 
theoretical limit which is rarely approached in practice. For instance, 
if one used this formula to calculate the separation resolvable by the 
human eye at 20 ft., it would be smaller than the actual minimum separation 
resolvable by a person with ^'normal" (20:20) eyesight. 

22 
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Problem 5 

The cones of the human eye happen to be maximally sensitive to the 
wavelength which corresponds to the maximum of the sun*s energy density 
spectrum. Using Rayleigh^s criterion and the fact that the retina is 
<r 17 mm behind the first nodal point (which acts as the circular aperture 
in this case), compute the maximum resolution on the surface of the retina. 

Answer : fix = 5.4 x lo"*^m. (13) 

Anatomists have determined that the diameter of a retinal cone cell 
(which mediates color vision) is about 3.0 x lO^^m. The separation 
.between densely packed cone cells is smaller than the, limit in Eq. (13). 
Thus thinner cones will not produce higher visual acuity. The resolution 
is constrained by the physical parameters of the optical system. 

Since the, eyes of all vertebrates are constructed from the same 
design, the same physical dimensions that determine the resolution of the 
human eye will determine resolution throughout the subphylum. It is the 
absolute, not relative, size of the physical components of the optical 
system which constrain resolution. If the diameters of cone cells are 
at present at the minimum which is biologically feasible., then a slight 
increase in the other dimensions of the human eye will give greater 
resolving power, but beyond a point (fix = 3.0 >^ 10 ^m), there is no 
payoff in visual acuity for larger eyes. 

The eyes of birds of prey show structural alterations which are an 
evolutionary response to the need for greater resolving power. The retinal 
cone cells become rod-like, which decreases their diameter, increases 
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their density on the retina and decreases their ability to discriminate 
color. Furthermore, the retinae from these birds have several pits or 
foveae. Although the density of light sensitive cells is the same in 
the foveae as on the rest of the retina, the image is formed on a concave 
rather than flat surface. Thus the "effective" cone density in the foveae 
is higher than for the rest of the retina, and so is the resolving power. 
The combination of these two changes results in increased visual acuity 
for these avian eyes. 

Rayleigh's criterion is also useful in elucidating the constraints 
on the design of a quite different sort of eye, the compound eye found 

in insects and certain Crustacea. Each eye is composed of many (up to 
4 

3 X 10 ) simple eyes, called ommatidia (See e.g. Bjorn 1976. p. 92). The 

shape of the compound eye is roughly that of a half-sphere of radius R, with 
the lenses of the ommatidia comprising its surface. An estimate of the 
size of the openings of the ommatidia is calculated in the following problem. 
Problem 6 

A simple model of an ommatidium is that of a cone with its apex at 
the center of the half -sphere which forms the eye, and its opening at the 
surface of the eye. Assume that the angle of the cone is equal to the limiting 
angle of the Raylcigh formula. Show that the calculated limiting size of the 
opening to the ommatidia is given as: 
d = (1.22RA)-^^^ 

If A = 400 nm, which is the center of the spectrum to which a bee is 
sensitive, and R is about 1 mm, what is the optimal size of the onimatidial 
opening? 

24 
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Answer: 

d « 2.2 X lO""^ m = 22 y. 

The actual mean measured opening for a bee is 20,9 y (Barlow 1952, p. 671) • 
A diagram comparing the observed size of the ommatidia to the Value predicted "by 
theory is presented in Fig. 4. 

Rayleigh's criterion for the limit of resolution of an optical instrument 
can be extended to interferometers as well. Basically, an interferometer 
consists of two detectors with the ability to detect phase infatlnation separated 
by a distance D, When the phase information from the two detectors is compared, 
the resolution of the total system is the same as that of a siivgle detector 
whose diameter is equal to D, Although a rigorous argument d^iftonstrating this 
result would not be appropriate here (it is essentially the saitte as the derivation 
of Rayleigh's criterion), a qualitative argument may be useful. Basically, the 
"'tiif fraction patterns illustrated in Fig. 3 are the result of tbe interference 
between light waves of different phases, caused by passage of the light through 
the aperture. The eye sees the interference caused by phase differences as 
variations in intensity, as illustrated in Pig, 3 b, c, d. In an Interferometer 
the signals of the two detectors are compared to reveal phase differences and 
so the identical resolution criterion applies. 

Most vertebrates are able to detect phase differences be»^veen ears, 
so that the interferometer is the correct model for calculating resolution. 
A consideration of Rayleigh's criterion shows that if all other parameters 
remain constant, the greater the separation of the detectors* the greater 
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FIGURE 4- INSECT EYE* OMMRTIDIRL DIRMETD^ VS. SORT. EYi: RADIUS 
POINTS' 27 OBSERVRTIONS ON HYMENOPTERfl (BARLOW 1952) 
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the resolving power of the system. Many organisms have optimized their 
ability to locate sound sources by maintaining a maximal separation 
between the ears. Another biological application* of interferometers is 
illustrated in the following problem. 

Problem 7 

Using the Rayleigh criterion and the. solution to Prob. 2b, determine 
the resolution of the infrared sensing system of a pit viper. Using 
reasonable physical parameters (e.g. body size of prey, how close the 
snake must be to strike, how accurate the strike has to be), determine 
the minimum separation of the pits. Does this agree with your intuition 
and/or experience? 

Answer ; d - separation of pits 6 cm. 

IV. The Doppler Effect 

The discussion up to this point has implicitly assumed that the 
source and detector of radiation are stationary relative to one another. 
When the source and detector are in relative motion, an added complication, 
the Doppler effect, is introduced. This phenomenon has been observed for 
all types of radiation, though the theory is simplest for the acoustic 
situation. Since acoustic waves are longitudinal (compression) rather 
than transverse as are electromagnetic waves, the mathematical form of the 
effect is slightly different. Additionally, the acoustic Doppler effect 
has significant evolutionary implications. 
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Empirically, the Doppler effect is manifestGd in the difference 
between the pitch (frequency) of the sound heard by a listener in motion 
relative to a sound source and the pitch heard when the listener and source 
are relatively stationary. A familiar example is the sudden drop in pitch 
one hears from an automobile horn as one meets and passes a car going in 
the opposite direction. Alternatively, one might think of the changing 
pitch heard by a listener standing beside a- track as a train whistle 
approaches > passes, and recedes. 

The derivation presented below is for the Doppler effect as it applies 
to the special case when the motion of the source and listener lies along 
the line joining the two. The biological example considered later is of 
precisely this type. 

Let and denote the velocities of the listener and source 
respectively. The positive direction for the velocities is taken to be from 
the position of the observer to the position of the source. 

An illustration of the situation when the listener L and the source 

S are moving away from each other is presented in Fig, 5. At time t = o, 

the source is at point X^. The outermost circle is the representation 

of the position of the wave front at time t = T, caused by a disturbance 

at the source at time t =. 0. The speed of propagation c of a wave in a 

s 

nondispersive medium such as air is dependent only upon the characteristics 

of the supporting medium and is independent of the motion of the source 

relative to the medium. Thus the outermost circle in Fig. 5 represents a 

sphere in 3 dimensions with center at x and radius C T. 

o s 

The source has moved a distance V T in the time T so that 

s 
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Figure 5: The Doppler Effect 




The source S, moving with velocity V , is emitting sound at a constant frequency, 

fg. The wave fronts of this^sound, which travel radially outward at velocity Cg, 

are compressed together (wavelength equals X ) in front of the source and spread 

F 

out (wavelength equals X ) behind the source. By the time the wavefront emitted 
by the source at Xq at time t = 0 has reached a radius of C^T, the source has moved 
to a position x^. The listner L moving with velocity will hear the sound at a 
different frequency than as actually emitted by the source because it encounters wave 
fronts which are X_ apart. An exact derivation is given in the text. [Adapted from 
Sears and Zenansky p. 327.] 
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and the following equalities are true: 

a - = (C + V )T 
T s s 



- b = (C - V )T 

T s s 

where a and b are respectively the positions at the rear and front of the 
outermost wave surface. 

In the time interval between t = o and t = T, the source has emitted 
a certain number of waves, f^T, where f^ is the frequency of the sound 
emitted at the source. The waves are spread out into the distance (x^ - b) 
in front of the source. 

In front of the source the wavelength is given by 
x^, - b (C - V )T 
s s 

Similarly behind the source the wavelength is 

a - (C + V )T 

X = ^ = — ^—-^ — = (C + V )/f . 

B f T f T ^ s s^ s 

s s 

The waves approaching the lijitener have a relative speed of + V^, 
so the detected frequency is 

C + 
s L 



S S 

or 

C + 

f = ^ ^ f . (1^0 

L C + V s 
s s 

When the source and listener are moving toward each other a frequency 

C + V, C + 

f = s L _ s L 



L X. (C - V )/f^ 

f s s s 
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or 



C + V 
f = s , 

L C - . (14a) 



Is detected. 

In order to help remember the direction convention established in 
this derivation, it is useful to notice that the frequency heard by the 
listener will be less than that emitted by the source if they are moving 
away from each other, and will be greater than that emitted by the source 
if they are moving toward each other. 

Example 5 

The velocity of sound in still dry air is = 350 m s'^ For a 
source emitting sound of frequency f^ = 700 E^ the wavelength of the 
sound emitted is 

^ = ^s'^'^s ° ^' 

a) What are the wavelengths of the sound in front of and behind 
this source moving at = 50 m;s~''"? 

b) If a listener is at rest and this source is moving away at 

= 50 n s what is the frequency of the sound heard? 

Splution: 

3) = ^^s ~ \^/^s =0-^29 m 
= ^°s \^/^s m 
C 

= C +% ■ = 612 Hz. 
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Now that a quantitative formulation of the Doppler effect has been 
derived, it can be applied in the explanation of echolocation in bats, a 
sensory system which has only begun to be understood in the past several 
decades. 

The bats, members of the order Chiroptera within the class Mammalia, 
have developed a system of echolocation which has been of immense 
evolutionary significance. It has enabled them to exploit a resource 
for which there are very few competitors, small night-flying insects 
(see, e.g., Fenton 1974). The' physical processes important in echo- 
location have both enhanced and constrained the evolution of the 
Chiroptera. Here it will suffice to show a few examples of the effect of 
- these physical processes upon bat evolution and natural history. For a 
detailed and personable account of the actual experimentation that led 
to the elucidation of the physical principles employed in echolocation, 
see Griffen*s Listening in the Dark . (Griff en 1958). 

Almost all bats appear to be able to navigate in total darkness a 
closely spaced grid of wires, where the spacings of the wires are 
commensurate with the dimensions of the wingspread of the bats. Furthermore, 
insectivorous bats can pursue and capture flying insects on the wing in 
total darkness. 

The ability to carry out these activities is significant in two respects. 
First, unlike the birds which are able to perform similar activities only in 
the presence of light, bats cannot be dependent upon a highly developed 
sense of sight for navigation. Second, no other animal can depend upon 
catching flying insects in the dark and so there are few competitors 
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for this food resource. Prior to the evolution of bats/ this resource 
had not been exploited. 

The experiments of Griffen demonstrate that bats share another 
ability- that sets them apart from most other animals. This is the ability 
to produce and detect sound of f requer. lies between 20-100. kHz» It is this 
attribute which is fundamental to the system of echolocation, as in the 
next example. 

Example 6 

A flying bat emits a pulse of sound. What is the minimum frequency 
and maximum length of a single frequency pulse, if the bat wants to use 
the echo of the pulse to avoid an obstacle 1 meter ahead? 

Solution: One first solves for the maximum pulse length. Clearly 
the chief consideration is that the bat must not be producing the pulse 
when the echo returns. Otherwise the fainter echo might be masked by the 
bat's own cry. The time of the echo's return t is dependent only on 
the speed of sound and the distance of the object d. Specifically: 

t = d/C 

s 

which works out to be: 

t = 2 m/350 m s"^ = 5.71 x lO"^ s = 5.71 ms. 

Shorter pulses will allow even closer objects to be avoided, although 
the pulse le.ngth must have a lower limit of several wavelengths in order 
to avoid problems in interpretation. If a reasonable number of waver, in 
the pulse is 100, the predicted minimum frequency would be: 
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f = 100/t = 100/4 X 10-^ s""-^ = 25 kHz. 



In actual field measurements Griff en (op cit% p. 191) has observed 



pulse lengths from 1-15 ms , and pulse frequencies between 30-75 k!: z 
with about 50 kHz being most comnon* Therefore, bats should be several 
times better at avoiding obstacles than one would expect from Example 6. 
This is not too surprising if the system is good enough to catch insects 
on the wing, which the bats must approach much closer than 1 meter. 

Interestingly, bats appear to be rather economical in their use of 
pulses, using relatively long pulses (15 ms) at long repetition 
intervals when they are cruising far d)ove the ground and relatively 
short pulses (1 ms) at short repetition intervals only when they are 
pursuing insect prey. 

Many other characteristics of the bat echolocation system can be 
explained by an extension of the kinds of arguments given here. For 
example, theoretically one expects that the intensity of the returning 
echo from an object at a certain distance, should decrease with the size 
of the object. The echolocation system provides information about the 
size of objects in the environment from the intensity of the returning 
echo, and infornation about the distance of objects from the time delay 
between the emission of the pulse and the reception of the echo. However, 
the character of the echo should theoretically be relatively insensitive 
to the detailed nature of the object causing the echo, whicli perhaps 
explains the interest exhibited by bats in pebbles tossed in the air by 
small boys and naturalists. 
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In certain bats, for example those of the genus Rhinolophus , the 
Doppler effect appears to be the physical phenomena which Is the key to 
an understanding of their system of echolocation I^Schuller, £t^ al. 1974), 
Several of the important considerations involved in the application of 
the Doppler effect to bats are developed in the following problem. 

Problem 3 

A bat is flying along at 5 m s It emits a 50 kHz sound pulse 
which is reflected by an insect which is 5 meters distant and is moving 
away at Im s What is the frequency of the echo, and how long after 
the pulse begins does the echo begin to return? 

Answer ; f = 51.74 kHz 

t = 28.6 ns 

The outgoing and Incoming frequencies i,n this example have a 
difference of 1740 Hz. This difference may appear quite large when 
compared to the difference between Concert A (440 Hz) and the A an 
octave higher (330 Hz). However, in the detection of radiation phenomena, 
the important comparison is the ratio, not the difference of two frequencies. 
The bat in Problem 8. must be able to detect 1740 Hz in 50 kHz, a percentage 
frequency change of 3.4%. This is not an unreasonable change to detect 
since human musicians can quite commonly distinguish quarter tones which 
represent a percentage frequency change of 3%. One expects that bats which 
depend upon the Doppler effect should have very acute pitch change detection 
and a very stable emission frequency, to go along with the ability to 
discriminate very small time intervals. Some of these predictions have 
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been experimentally verified (Schuller, eit al . op cit Q, 

The thrust of the discussion presented in this section is perhaps 
most striking when viewed in terms of the advantages and constraints 
that the utilization of echolocation has conferred on the evolution of 
the Chiroptera. The development of echolocation has allowed the bat to 
utilize a previously unexploited resource, night flying insects, but in 
turn has placed a premium on those bats that are able to emit high 
frequencies in short pulses, that caA detect faint echoes, and that can 
detect relatively small shifts in frequency between outgoing and. incoming 
pulses. The selection pressure for these attributes is the force respons- 
ible for numerous remarkable physical and behavioral adaptations. Our 
understanding of this system is by no means complete; echolocation remains 
a fertile territory for experiment and theory. 

V. SUMI-IARY 

The preceding discussion has presented some of the results from 
classical and modern radiation theory and demonstrated how these physical 
principles can be applied to the analysis of biological systems. The 
emphasis in these applications has been upon the evolutionary significance 
of the physical processes. In the biological applications, presented, the 
physical processes are involved in the following two types of explanations 

In the first type of explanation; physical theory is used to isolate 
constant factors in the environment to which all organisms must adjust 
their evolution. The most striking example of this is that the wavelength 
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RADIATION 

Appendix 1. Symbols, Units and Dimensions 



Symbol Quantity 



Units 



Dimension 



a 
b 
c 
C 

3 

D 
d 

^i 

dX 

dv 
56 
fix 
^ 



E(X) 
E(V) 



angle radians 
Wien constant 
speed of light 
speed of sound 
distance 

distance to an object 
i^^ difference term — 
differential wavelength m 
differential frequency 
infinitesimal angle radians 
infinitesimal distance m 
total power J s 



power flow at 
earth's surface 

radiation emittance 

power flow at 
sun's surface 



2.8978 X io""\ X 
2.998 X lo^m s""^ 
3.50 X 10^ m s^"^ 



m 



m 



-1 



-1 



T -2 -1 

J m s 

T -2 -1 

J m s 

T -3 -1 

J m s 



H 



energy density as a J m 

function of wavelength 

energy density as a 
function of frequency 

base of natural logaritlims 

frequency at listener 

frequency at source 

total heat loss rate J 5 



-4 



J m"^ s 



2.171828... 
-1 

s 

-1 

s 

-1 



LS 

lt' 

LT 
L 
L 



-1 



rl 



2 -3 
ML T 




-2 -2 



-1 



-1 



2 -3 
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Appendix I (cont'd) 



Symbol Quantity 



Units 



Dimension 



Hz 
h 
k 
X 



0 
R. 



s 

s 

o 
SA 



Hertz 

Planck constant 
Boltzmanri constant 
wavelength 

wavelength behind source 

wavelength in front of 
source 

wavelength at which black- 
body spectrum has maximum 
energy density 

observer 

i^^ remainder term 
radius 

mean distance from earth 
to sun 

radius of sun 
source 

ft 

size of an object 
surface area 

Stefan-Boltzmann constant 

temperature 
a particular time 
ambient temperature 
temperature of organism 
surface temperature of sun 
time variable 



-1 



6.626 X 10" s 

-23 -1 
1.381 X 10 ^ J K 



m 

m 
m 

m 



m 



7.1 X io*m 



m 



-8 



5.6696 X 10 

-1 -2 -4 

J s m K 



K 

s 
K 
K 
K 



-1 



2 -2 -1 

ML T e ^ 



L 
L 
L 



L 
L 



T 

e 
e 

9 

T 
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Appendix 1 (contM) 



Symbol Quantity Units Dimension 



L = length 



M = mass 
T « time 



6 = temperature 



u actual solution 

V frequency s T 

V- listener velocity m s T ''"L 

Li 

source velocity m s 

i^^ trial solution — — 

X * point in space corres- m 
° ponding to source position 

at t = 0 

x^ point in space corres- in 

ponding to source 
position at t = T 
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Appendix 2: I. Numerical Techniques 

Introduction 

In the discussion of black-body radiation there arose two equations 
which could not be solved analytically. Since such equations arise not 
infrequently in a variety of areas, for example in mathematical modeling 
of biological systems, it is worthwhile to investigate methods by which 
one can find solutions to such equations. One class of such methods, called 
numerical techniques, utilizes the recursive application of the basic oper- 
ations defined on the real number system (+, -,. to find solutions 
accurate to any desired accuracy. With the increased availability of 
machines which can carry out these techniques swiftly and automatically, 
the problem solving ability of the scientist has Increased enormously. In 
this appendix, two important numerical techniques, successive approximation, 
and the half increment method for definite integrals, will be presented and 
their implementation on the increasingly ubiquitous, programmable pocket, 
calculator will be demonstrated. 

Successive Approximation 
As an example of this technique, the solution of the transcendental 
equation 

e"" + ~u - 1 = 0 (A-1) 
which arose in the discussion of Wien's Law will be discussed. 

Here u is the unknown to be calculated. Suppose one starts with a 
trial solution x^ . If d^ is the difference between the first trial 
solution X, and the actual solution u, then 
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X, = u - d^. 



Substituting for u in eq. (A-1) gives: 

-X; 

5^^! - ~ 



e ^ + -Ix - 1 = R, 



a-2) 



th 



where Rj is the first remainder term. Clearly as R^ -> 0 the n trial 

solution will approach u. 

Subtracting eq. (A-1) from eq. (A-2) one obtains: 

"X, -u. 



^1 -1 
e - e 



or 



-X, 



-d. 



1 - e 



(A-3) 



The exponential in the unknown dj in eq. (A-3) can fee expanded as an 
infinite series (e.g., a Taylor series): 



-d, d,^ d,-^ 



If the first trial solution is a reasonable guess then [dj | will be 
small so that the first few terms in the infinite series will dominate. 
Using only the first two terms, eq. (A-3) becomes: 



""^1 1 



or 



dj (e - 1/5) = Rj 



or finally: 



= R /(e ' - 1/5) 
Clearly the second trial solution x^ , defined by 



(A-4) 



Substituting for u in eq. (A-1) gives: 



e + -^x^ - 1 « a-2) 

where Rj is the first remainder term. Clearly as R^ -> 0 the n^^ trial 
solution will approach u. 

Subtracting eq. (A-1) from eq. (A-2) one obtains: 



"Xj -Uj 

e - e 



-Id, =R, 



or 



- J dj = Rj (A-3) 



The exponential in the unknown dj in eq, (A-3) can fee expanded as an 
Infinite series (e.g., a Taylor series): 

-d, d,^ d,-^ 

e =i-di+-2r--3r"*'Tr" 



If the first trial solution is a reasonable guess then [dj | will be 
small so that the first few terms in the infinite series will dominate. 
Using only the first two terms, eq. (A-3) becomes: 

or 

dj (e - 1/5) = 

or finally: 

dj = R^ / (e - 1/5) (A-4) 
Clearly the second trial solution , defined by 
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X, = u - d^. 



Substituting for u in eq. (A-1) gives: 

-X; 

5^^! - ~ 



e ^ + -Ix - 1 = R, 



a-2) 



th 



where Rj is the first remainder term. Clearly as R^ -> 0 the n trial 

solution will approach u. 

Subtracting eq. (A-1) from eq. (A-2) one obtains: 

"X, -u. 



^1 -1 
e - e 



or 



-X, 



-d. 



1 - e 



(A-3) 



The exponential in the unknown dj in eq. (A-3) can fee expanded as an 
infinite series (e.g., a Taylor series): 



-d, d,^ d,-^ 



If the first trial solution is a reasonable guess then [dj | will be 
small so that the first few terms in the infinite series will dominate. 
Using only the first two terms, eq. (A-3) becomes: 



""^1 1 



or 



dj (e - 1/5) = Rj 



or finally: 



= R /(e ' - 1/5) 
Clearly the second trial solution x^ , defined by 



(A-4) 
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TABLE A-1 (SR-52) 



Solution to e~" + -ju - 1 



00 



01 

X. 



02 



Register Contents: 
03-18 19 



Used for calculation of R^SBR 1*) 



Program; 



Step 


Code 


Key Entry 


Comments 


000 


42 


STO 


in display 


001 


00 


0 




002 


01 


1 


stored In register 01 


003 


51 


SBR 




004 


87 


1* 


Calculates 


005 


42 


STO 




006 


GO 


0 




007 


02 


2 


stored in register 02 


1)08 


81 


HLT 


Cease execution, display 


009 


43 


RCL 




010 


00 


0 




Oil 


01 


1 


Recall 


012 


51 


SBR 




013 


23 




-X. 

Calculates e ^ 


014 


75 






015 


93 


• 




016 


02 


2 




017 
018 


95 
20 


l/x 


"\ 1 

e - in display 

"^i 1 

(e - in display 


019 


65 


X 




(Continued) 









Note 1** refers to 2nd 1. 



47 



TABLE A-1 (SR-52), Continued 



Step 


Lode 








020 


43 








021 






u 


• 


022 


02 




Z 


i\eCa±± Iv^ 


023 








r1 St P ^o~^i -fn rl-fcnliJV 


024 


oc 

85 




t 

T 




025 


43 




DPT 




026 


00 




u 




027 


01 




1 

1 


Kecaii 


028 


oc 

95 






X. , , ^ X, T a xii uxc»pj.ay 
i+1 1 , i 


029 


81 




HLT 


Cease execution, display ^^^.^ 


030 


OD 




1. 0 O U 


pot-fiT*r» t-o et*oT\ nnn 


•^31 


4o 




T "RT 




032 


23 




0 rw 


DUDX^OULHie XllX 


033 


94 




+/- 




034 










035 


oo 




XrllX 




036 


56 




rtn 




03/ 


/.A 
4o 




T VI 
LiIjLi 




038 


87 




1* 


Sul;* -tine 1* 


039 


53 




( 




OAO 


42 




STO 




041 


01 




1 




042 


09 




9 


Store in register 19 


(Continued) 




NOTE: 


1* refers to 
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TABLE A-1 (SR-52), Continued 



Step 


Code 


Key Entry 


Comments 


043 


51 


SBR 


• 


044 


23 


Jtnx 


Calculate e"^^ 


045 


75 


— 




046 


01 


1 


Calculate e 


047 


85 


+ 




048 


53 


( 




049 


43 


RCL 




050 


01 


1 




051 


09 


9 


Recall 


052 


55 


• 




053 


05 


5 


Calculate ^ 


054 


54 


) 




■t55 


54 


) 


Calculate 


056 


56 


rtn 





i9 
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B. The HP-25 

Using the Program (Table A-2): 

Switch the PGRM-RITN switch to PGRM. Clear by keying 
f PGRM . Then key steps 01 through 24 of the program. Switch PGRM-RUN 
switch to RUN , and set the calculator to step 00 by keying f PGRM . Key 
Xj into the display. Run the program by pushing the R/S button. The 
calculator will cease executing when R^^has been calculated. R^ will be 
in the (x) display. Push R/S ; after calculation x^, will be displayed. 
Push R/S; the display will show R^, etc. 

To try a different first trial solution, key f PGRM and then 
key in the new x^. Run the program as before. 

There are two solutions to eq. (A-1): u = 0.0 and u - 4.9651... 

If Xj is too small or negative, this program will converge to the trivial 

solution (u = 0.0). Try a larger Xj ; i.e., key f PGRM , and then key in 

the new x^. Run the program as before. If x^ is too large, the display 

will show OF when the calculator ceases execution. Start over again with 

a smaller x . 

1 

For troubleshooting and editing of a HP-25 program, see section 
IV .B. of this appendix. 
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(Continued) 
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TABLE A-2 (HP-25)'. , Solution to Transcendental Equation 

Register Contents 
0-7 0 



Step 


Code 




X 




Z 


T 


00 












01 


23 01 


STOl 








■ 


02 


31 


ENTER 


h 


h 






03 


32 


CHS 




h 






04 


15 42 


X 

ge 










05 


31 


ENTER 




-x, 

e 






06 


31 


ENTER ; 


e * 


e 




Xj 


07 


73 


t 


* 


-Xi 

e ^ 


-X, 

e 1 


Xj 


08 


02 


2 


.2 




e 1 


^1 


09 


32 


CHS 


-.2 


-Xi 

e * 


e 1 


Xj 


10 


51 


+ 


e"''-.2 


e * 


h 




11 


23 02 


ST02 


e 1-.2 


e I 






12 


22 




e-^i 




\ 


e l-,2 


13 


21 


x<y 


h 









Copent 
Xj in display 
Start loop 
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TABLE A-2 (HP-25). Continued 



Ste£ 


Code 


Key 


X 


JL 


Z 


T 


Conroent 


14 


31 


ENTER 






e-xi 


h 




15 


05 


5 


5 


^1 




h 




16 


71 


t 


1/5 xi 


"Xi 

e 1 




h 




17 


31 


ENTER 


1/5 Xj 


1/5 Xj 


e 1 


h 




18 


01 


1 


1 


1/5 xi 


"Xi 

e * 


h 




19 


41 




1-1/5 


-X 

e 1 




h 




20 


51 


+ 




Xj 




h 


Calculate \ 


21 


74 


R/S 


h 


Xj 


h 


h 


Display Rj 


22 


24 02 


RCL2 




Rj 


h 


h 




23 


71 


J. 
• 


di 


X] 


h 




Calculate dj 


24 


24 01 


RCLl 


^1 


dj 


h 


h 




25 


51 . 


+ 






h 


h 


Calculate Xj 


26 


74 


R/S 


Xj 


x^ 


h 




Display Xj 


27 


13 01 


GTO 01 


Xj 


Xj 


h 


h 





^ 

V 

o 

ERIC 
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C. Using Successive Approximation to Solve Other Problems 

The method of successive approximation is a powerful technique for 
finding the solution to certain equations quickly and easily, as illustrated 
by the programmed solution to eq. (A-1) . Besides transcendental equations 
(viz. (A-l)), this technique is useful for solving: 

1) Poljmomials with constant coefficients of any degree 

2) Polynomials with variable coefficients of small degree. 
The crux of the method ^3 finding an explicit formula for a 

correction term (d,^^) so that from any trial solution (x^^^) one can calculate 
a new solution (x^^^) which is a better approximation to the correct 
solution than the original (x^). To find this correction formula often 
strains the ingenuity of the problem solver, who may have to juggle 
equations or approximate functions by series. A particularly useful method 
for finding these correction formulae, especially for the polynomial 
problems listed above, is the Newton-RaptiSon method which is derived and 
discussed in almost all introductory books on numerical methods^. 

Definite Integrals by the Half- Increment Method 

In the derivation of the Stef an-Boltzmann law, it was necessary to 
evaluate the indefinite integral: 

CO 

-7^ dx. (A-5) 
•'o e -1 

Hsee, e.g., Weeg & Reed 1966, p. 34.] 
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This seemingly simple integral cannot be evaluated analytically i Unfortu- 
nately, this is not an isolated example. That is, the majority of functions 
including many important simple functions, cannot be integrated by standard 
analytical techniques. In such cases, one must resort to numerical m^ithods 
to obtain the desired integral. The particular numerical integration 
technique described in this appendix, the half-increment method, is simple 
to understand and to implement on a programmable pocket calculator, and has 
the added advantage of (relatively) rapid convergence. 

In elementary calculus, one learns that an integral can be thought of 
as a measure of an area. For a definite integral such as: 

b 

I f(x)dx (A-6) 

a 

this area is defined by the curve of the function to be integrated, f(x), 
the X-axis, and the vertical lines x - a and x = b. This area can be 
approximated by the sum of the area of rectangles as illustrated in 
Figure A-1. The thinner the rectangles and thus the greater the number 
of rectangles, the more accurate the approximation. In fact, it can be 
shown that in the limit as the width of the rectangles becomes infini- 
tesimal (Ax 0 in Fig. A-1) , the value of the sum approaches the value 
of the integral. 

All numerical integration techniques are based on approximating an 
integral by the sum of the areas of simple plane figures. This particular 
method is called the half-increment method because the height of each 
rectangle is determined by the value of the function at the x value which 
is half-way between the two x values which define the vertical boundaries 
of the rectangle. That is, for the i*''^ ^rectangle, whose vertical 

56 
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Figure Definite integral approximated by rectangles 




The value of the integral, that is the area under the curve f (x) from x = a 
to X = b, is approximated by the sum of the area of the rectangles illustrated 
above. In this approximation, the height of each rectangle is taken to be the 
value of the function at the x value halfway from the endpoints of the rectangle 
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boundaries are determined by x = and x = = x^ + Ax, the height of 

the rectangle is given by f(xj^ + = f (x^^^^^ - • 

In the programs that follow, the integral (A-6) is evaluated first as a 
sum of rectangles whose width, A^x, is given by A^x = b - a. This sum, 
Ej, which is composed of only one term, is given by: 

A,x 



AjX 



Next the sum, E^, is calculated for A^x = -y . At each succeeding 
iteration, the width of the rectangles is halved; that is: 



and 



2 

3-1 

where 



i f 



A.x 

1 



= a and x ^ = b. 

2 

When the value of the sum as determined by two successive iterations is 
the same, this sum is taken to be the value of the integral. Obviously, 
each iteration will require twice as many operations as the preceding, 
and will take correspondingly longer. 

The integral of interest, (A-5), is indefinite, although the method as 
presented is only applicable to definite integrals. Fortunately, the sum 
does not have to be evaluated for very large b (the upper limit in (A-6)) 
in order to be accurate to two decimal places. This is a result of the 
following limits: 

58 
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3 3 

(which can be calculated by L'Hospital's rule), and the fact that the 
function approaches zero very quickly (and monotonically) for sufficiently 
large x. 

A. The SR-52 

Using the program (Table A-3): 

With the calculator in the execute mode, key LRN to enter program- 
ming mode. Key in steps 000-111. Key LRN to return to the execute mode. 
Push rset to go to step 000. Load the value of a in register 01 (i.e., 
key in the value of a, aad push STO 01,), and the value of b in register b. 
Push RUN . Vl.an the calculator ceases execution of the program, will be 
displayed. Push RUN. 2^ will be displayed, and 2^ will be in register 06. 
Run the program until 2^ "^i-i ^-^i ^^^^ in register 05, l^^^ will be 
In register 06). Note that each iteration takes twice as long as the 
preceding one. 

Clearly the lower limit of the sum should be taken as the same 
as the lower limit of the integral (A-5); that is, a *= 0.0. However, the 
upper limit (b) is not as easy to choose. The larger the value of^ b, the . 
better the approximation of the sum to the int<5gral.. However, concomitant 
with a larger b is a greater number of iterations required for convergence. 
A good compromise value is b = 30, which will approximate the actual value 
of the integral (I = 6.4938...) to 3 decimal places in 9 iterations. The 
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program user should try other (larger and smaller) values of b In order 
to get some feeling for these tradeoffs. 

For troubleshooting and editing of an SR-52 program, see 
Section IV. A of this appendix. 



/ 
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TABLE A-3 (SR-52). Integral by Half-increment Method 
Register 



00 



01 
a 



02 
b 



preloaded 



03 



04 
Ax 



05 



Contents 
06 
^i-1 



07 - 18 



19 

used for f(x) 



Program 



Step 


Code 


Key Entry 


Comments 


000 


00 


0 




001 


42 


STO 




002 


00 


0 




003 


05 


5 


Zero register 05 


004 


43 


RCL 




005 


00 


0 




006 


02 


2 


Recall b 


^07 


75 






008 


43 


RCL 




009 


00 


0 




010 


01 


1 


Recall a 


Oil 


42 


STO 




012 


00 


0 




013 


03 


3 


Store = a in Register 


014 


95 




Calculate A^x = b-a 


015 


42 


STO 




016 


00 


0 




017 


04 


4 


Store A-x in Register 04 



(Continued) 
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TABLE A-3 


(SR-52), 


Continued 




Step 


Code 


Key Entry 


Comments 


018 


55 


* 


• 


019 


02 


2 




020 


95 


- 


Ax 

Calculate 


021 


85 


+ 




022 


43 


RCL 




023 


00 


0 




024 


03 


3 


Recall X. 

1 


025 


95 




Calculate x. + -5- 


026 


51 


SBR 




027 


00 


0 


This subroutine 


028 


09 


9 


calculates 


029 


05 


5 


f (x^ f^) 


O30 


65 


X 




031 


43 


RCL 




032 


00 


0 




033 


04 


A 


Recall 


034 


95 


S3 


Ax 

Calculate f (x^ + -y) 


035 


44 


SUM 


Put 


036 


00 


0 


^Ef (x^ + ^) Ax in 
i 2 










037 


05 


5 


Register 05 


038 


43 


RCL 




039 


00 


0 





S2 

(Continued) 

o 

ERIC 



TABLE A-3 CSR-52), Continued 







xvey £#n t ry 


Coimuents 


\JH\J 


HA 




Recall Ax 




AA 


SUM 


uaxcuxace 


HA 9 


uu 


n 
u 


''i+i = + ^ 




n't 




Store in Register U3 




A*i 


T>r»T 
KCL 




HA *i 


uu 


U 












ftA 7 


7Q 
/3 






HAQ 


AQ 


KLXj 




HAQ 


uu 


u 




f\Kn 


no 
u^ 


9 


TJ #» 1 1 "k 

Kecaxx D 








L>£i±cuxaL.e ^ji^^j^ D 




fin 

oU 


XI. jL/UO 


i-fl ^ ^ SO 


UjJ 


^o 


SXIl 


seep xaDe±ea ^ 


n*^A 


A^ 


KLXj 


If ^i+l < b 


0!>!> 


nn 
UU 


U 


Kecaxx 


056 


04 


4 


Ax 


057 


41 


GTO 


Continue 


058 


00 


0 


the sum 


059 


01 


1 




060 


08 


8 




(Continued) 
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TABLE A-3 (SR-52), Continued 





Code 


Kev Entrv 


061 


46 


LBL 


062 


32 


^ "i n 

O Jul! 


063 


43 


RCL 

XV W li< 


064 


00 


0 


065 


05 


5 


066 






067 


43 




06ft 


00 


0 


06Q 


05 


5 


070 


42 


STO 


071 

U/ J. 


00 


0 








079 


06 

UU 


6 


073 


00 


0 


07A 


42 


STO 


07 S 


00 


0 


076 


05 


5 


077 


02 


2 


078 


22 


INV 


079 


49 


PROD 


080 


00 


,0 


OfU 


04 


4 


082 


43 


RCL 


0S3 


00 


0 



Comments 



Recall \ 
Display \ 



Recall I 



Store \ in Register 06 



Zero Register 05 (I) 



A.x 

Put Aji^^-j^'^ "T~ 

Register 04 



continued 
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lABLE A-3 (SR-52), Continued 



continued 



Step 


Code 


Key Entry 


084 


01 


1 


085 


42 


STO 


086 


00 


0 


087 


03 


3 


088 


43 


RCL 


089 


00 


0 


090 


04 


4 


091 


41 


GTO 


092 


00 


0 


093 


01 


1 


094 


08 


8 


095 


42 


STO 


096 


01 


1 


097 


09 


9 


098 


45 




099 


03 


3 


100 


55 




101 


53 


( 


102 


43 


RCL 


103 


01 


1 


104 


09 


9 


105 


22 


INV 


106 


23 


Anx 



S5 



Comments 



Recall a 



Store = a in Register 03 



Recall A , , - X 
i+1 

Return to Step 018 
to recalculate 

Integral w/ Smaller 

Ax. 

Store 

, Ax . 

Register 19 



Calculate 



r ^ Ax] 

X. + -TT 



Recall 



Calculate exp 
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TABLE A-3 (SR-52), Continued 



Step 
107 
108 
109 

110 
111 



Code 
75 
01 
54 

95 
56 



Key Entry 

1 
) 



Comments 



Calculate exp|x^ + -^j - 1 



Calculate f 



rtn 
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ERIC 
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B* The ET 

Using the program (Table A-4): 

Switch the PGRM-RUN switch to PGRM, Clear program memory by pushing 
f PGRM , Then key in steps 01 through 46 of the program. Switch PGRM-RUN 
switch to RUN, and set the calculator to step 00 by keying fJPGRM. Load a 
in register 1 (i.e., key a into the display then push STOJL) and b in 
register 1. Push R/S > The display will show Z^, Push R/S. The display 
will show Sj. will be in register 6. Run the program until (in the 
display at the end of execution and in register 5) equals ^ (in register 
Note that each, iteration takes about twice as long as the previous one. 

Clearly the lower limit of the sum should be the same as the lower 
limit of the integral (A-S); i.e., a = 0.0. However, the upper limit of 
the sum (b) is not as easy to choose. The larger the value of b, the 
better the approxiroiation of- the sum to the integral, but the longer it 
takes for two successive iterations to converge. A good compromise value 
is b = 30, which will approximate the actual value of the: integral 
(I = 6.4938.*.) to three decimal places in nine iterations. The program 
user should try other (larger and smaller) values of b in order to get 
some feeling for the tradeoff between accuracy and time in using this 
method on this integral. 

For troubleshooting and editing of an HP-25 program see section IV. B 
of this appendix. 



67 



TABUS A-4 (HP-25) 



a 



■ — dx by Half-Increment Method 
e -1 




Register Contents 
4 5 
Ax I 



i-l 



preloaded 



Program! 



Step 


Code . 


Key Entry 


X 


Y 


01 


00 


0 


0 




no 


23 05 


STO 5 


0 




03 


24 02 


RCL2 


b 




04 


'24 01 • 


RCLl 


a 


b 


05 


23 03 


STO 3 


a 


b 


06 


41 




AjXab-a 




07 


23 04 


STO 4 


Ax 




08 


02 


2 


2 


Ax 


09 


71 


T 


Ax/2 




10 


24 03 


RCL3 


^i 


Ax/2 


11 

(Continued) 


51 


+ 


x^+Ax/2 


V 



Conments 



Zero register 5 
Recall b 
Recall a 



Calculate A^x»b-a 



Continue loop [from step 36] 
Calculate Ax/2 



Calculate x^+Ax/2 



"TO" 
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BLE A-A (HP-?5), Oontinutid 



2A 04 
61 

23 51 05 
2A 04 
24 03 
ontlnued) 

\ 



RCL 4 

X 

STO+5 
RCL 4 
RCL 3 



Ax 

f(x^+Ax/2)A:: 
f(x^+Ax/2)Ax 
Ax 



f(x^+Ax/2) 



f(x^+Ax/2) 



Ax f(x^+Ax/2) 



Code 


Key Entry 


_X 


Y 


Z 


T 


Comments 


31 


ENTER 


x^+Ax/2 


x^+Ax/2 






Calculate exp[x^+Ax/2] 


15 07 


X 


exp[x^+Ax/2] 


x^+Ax/2 








01 


1 


1 


exp[x^+Ax/2] , 


x>Ax/2 






^1 




exp[]-l 


x^+Ax/2 






Calculate exp[x^+Ax/2] 


•'.1 


> 

X < y 


x^+Ax/2 


exp[]-l 








31 


ENTER 


x^+Ax/2 


x^+Ax/2 


exp[]-l 






03 


3 


3 


x^+Ax/2 


exp[]-l 






14 03 




(x^+Ax/2r 


exp[]-l 






Calculate (x^+Ax/2)^ 


21 


> 

x < y 


exp[]-l 


(X.+AX/2)' 






(:;.+Ax/2)5 


71 




f(x^+Ax/2) 








f(x^+Ax/2) = 



exp[x^+Ax/2]-l 



Calculate f (x^+Ax/2)Ax 
Sum cf areas in reg. 5 
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Table A-4 (IIP-25), Continued 



Ste£ 


Code Key Entry 




JL 


J. 


-L 


Comments 


27 


51 + 


1 


f(x^+Ax/2) 








28 


23 03 ST0 3 


^+1 


f(x^+Ax/2) 




• 


x^^^ in reg. 3 


29 


24 02 RCL2 


b ■ 


''i+l 


f(x^+Ax/2) 




Test for x , beyond 
l+i 


30 


14 41 fx<y 


b 


^+1 


f(x^+Ax/2) 




range of integration ' 


31 


13 36 GT0 36 


b 


^+1 


f(x^+Ax/2) 




of Xj^^j^>b, leave loop, 



32 24 04 RCL4 

33 31 ENTER 

34 13 08 GT0 08 

35 24 05 RCL5 

36 74 R|s . 

37 23 06 ST0 6 

38 24 01 RCLl 

39 23 03 ST0 3 

40 • 14 34 f STK 

41 25 05 ST0 5 



Ax 
Ax 
Ax 
I 
I 
I 
a 
a 
0 
0 



b 


^i+1 


f(x^+Ax/2) 


Ax 


b 


Vl 


Ax 


b 


^i+1 


b 


^i+1 


f(x^+Ax/2) 


b 


Vl 


f(x^+Ax/2) 


b 


^i+1 


f(x^+Ax/2) 


I 


b 


*i+l 


I 


b 


Vl 


0 


0 


0 



go to step 36 



If < b 



Continue loop at step 08 



Cease execution, display I 
l^^l^,^ in reg. 6 

Xj^ B a in reg. 3 
Zero stack 
Zero reg. 5 



(continued) 



EMC ' 
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TABLEAU (HP-25), Continued 



Step 


Code 


Key Entry 


X 


Y 


42 


02 


2 


a 




Hi 


91 71 n/i 

a IL Uh 




i 






24 04 


RCL4 


Ax 


2 


45 


il- 


X ^ y 


2 


Ax 


46 


ia 09 


GT0 09 


2 


Ax 



Coiments 



lower limit from reg 2 
Ax- 



Recalculate sum with 
smaller rectangles 



vD 
\0 



74 



75 



ERIC 
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C. Using the Half- Increment Method to Integrate Other Functions^, 
The SR-52: 

The subroutine 095 calculates the value of f(x^ integrate 
any other function^ define the fxmction with suitable keystrokes starting 
at program step 095 sure to end the fiinction definition subroutine 

with rtn. 

If the program is already loaded into program memory it is possible 
to change only those statements involved in the functional definition 
(subroutine 095), which is usually less trouble than wiping program memory 
and starting anew. For appropriate editing instructions, see section II. A 
of this appendix. 

The HP-25: 

The steps 12 through 22 are used to calculate the value of f (x^ + . 
To integrate any other function, define the fimction with suitable key strokes 

Starting at step 12. Be sure to end the sequence with f (x^ + ^) in the 

X register. Follow the last step used to define the function with the step 

numbered 23 in the pr^^sent program (that is, the keystrokes RCL 4) . 

If the program is already loaded into program memory, it is possible 

to change only thO£> keystrokes involved in the functional definition 

(steps 12-22) . This is usually less trouble than wiping program memory and 

starting anew. For appropriate editing instructions see Section II. B of 

this appendix. 



^ [For alternative programs and more detailed discussion, see 
ErJc Elsberg, R. M. 1976.3 

?6 
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Appendix 2: II. Progran Troubleshooting and Editing 

It is frequently the case that after a program has been keyed into a 
calculator, the program will not perform as expected. This frustrating 
state of affairs may be due to a logical error in the program or mistakes 
made when keying in the program; e.g., leaving out steps, putting in extra 
steps, or keying steps in the incorrect order. Since such mistakes are 
relatively common, the design of programmable pocket calculators is such 
that troubleshooting and editing programs is not difficult. 

Mistakes involving logical errors, which are the most difficult to 
detect, will not be present in the programs presented here. For such errors 
in your own programs, I have found that "single-stepping" through the programmed 
calculations (as described below) is the most efficient method of detecting 
the error. For further information, refer to the owner's handbook supplied 
with the calculator. 

For nonlogical programming errors, there is always the option of wiping 
the program memory (usually by turning the machine off) and starting over. 
This procedure, besides being time consuming^ does not assure that additional 
errors will not be made in loading the program the second time. Alternatively, 
one can edit the program as described ijelow. In the following discussion, it 
is assumed that the program is already loaded into program memory. 

A. The SR-52 

Troubleshooting: 

Programming errors are indicated by any of three conditions following 
a program RUN : 
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1) The calculator does not cease execution In a "reasonable" length 
of time. This condition Is due to a closed loop In the program instructions. 
The calculator cannot exit the loop so as to reach a HLT instruction. 

2) Vhen the program ceases executing, the display is steady, but the 
answer is clearly incorrect. This situation ©ay be difficult to detect 
unless one knows exactly what the program does. 

3) When the program ceases executing, the display is flashing. This 
condition signifies that during execution, a calculation was performed that 
overflowed the capacity of the calculator, or that an illegal operation was 
performed. 

There are two methods by which to pinpoint the source of an error. One 
can check the program in the memory against the program 'listing, or one can 
single-step manually through the calculations that the calculator performs 
automatically. 

To check the program in the memory, push rset. and then LRN to enter 
programming mode. The display will show the key code and step number of 
step 000. Pushing SST will display the key code and step number of the next 
step. Pushing bst will display this information for the preceding step. 
Compare the key codes displayed to those given in the program listing. If 
there is some reason to believe that the error occurs at some particular 
step, e.g.. Step 099, one can go immediately to this step, in the execute 
mode by pushing GTO 099 , Then push LRN to enter the programming mode and 
to display the key code and step number 099. Use SST and bst as necessary. 



Alternatively, one can single step through the programmed calculation 
manually. This is particularly useful for debugging your own programs or when 
a flashing display is encountered.. In the execute mode, push rset S^t 
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to step 000. Load the registers and display as necessary. Continue in the 
execute mode and push SST. The calculator will execute instruction 000 and 
move on to the next step. Pushing SST repeatedly, one can observe the result 
of each programmed instruction as it is performed. If at any step an error 
condition results, push LRN to enter programming mode, and then push bst. 
The step number and key code of the incorrect step will be displayed. 

Editing: 

Deleting instructions: In the program mode, with the step number and 
key code of the incorrect step displayed, press del. The incorrect instruction 
will be deleted and the next instruction will be shifted down to this step 
number, and all subsequent instructions will shift down one step number. 

Changing instructions: In the program mode, with the step number and 
key code of the incorrect step displayed, press the correct key for that 
stepnumber. The correct instruction will write over the previous one, and 
the display will shift to the next step. 

Adding instructions: In the program mode with the step Tiumber and key 
code of the step displayed, press INS. The display will change to the correct 
step number but the key code will be 00. The key code that was at this step 
will be shifted to the next step number, as will all subsequent instructions. 
Press the key for the missing instructions. The display will now show the step 
number and key code of the next step. 

B. The HP-25-. 

Troubleshooting : 

Progra mm i n g errors are indicated by any of four ct. ditions following a 
program run: 
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1) The calculator does not cease execution after a "reasonable" 
length of time, this condition is due to a closed loop in the program 
instructions. The calculator cannot exit the loop so as to reach a 

R| S instruction. 

2) When the program ceases executing, the answer displayed is 
clearly incorrect. This situation is difficult to detect. 

3) When the program ceases executing ERROR is displayed. This 
results from attempting an illegal instruction. The program stops at the 
incorrect instruction, which will be displayed by switching to PGRM code. 

4) When the program ceases executing, OF is displayed. This results 
from an overflow in a storage register. Switch to PGRM mode to display the 
key code and step number of the instruction that caused this condition. 

To compare the program in memory to the program listing, in the RUN 
mode push GTOOO (or if there is reason to believe the error is at some 
later step, e.g., step 16, push GTO 16 . Pushing SST will display the step 
number and key code of this step. When SST is released, the instruction 
is executed, the result is displayed, and the calculator moves on to the 
next step. 

Editing: 

Changing an instruction. Once the error is located, use the GTO (in^ 
RUN mode), the SST, or the bst (in PGRM mode) necessary so as to be at 
the step preceding the step which is to be changed. In the PGRM mode, with 
the key and step number of the preceding step displayed push the key of 

the correct instruction. 

Deleting an instruction: Display the step preceding the step to be 
deleted (using GTO , SST, or bst as necessary) . In the PGRM mode press 
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8 NOP , When the calculator reaches this program step, it will skip it and 
proceed to the next instruction. 

Adding instructions: To add instructions, one uses the GTO instruction 
to perform an unconditional branching* For example, suppose one has left out 
three steps between Step 08 and 09. Further, suppose the program ends with 
Step 15. Then one replaces Step 08 with GTO 16 . Then key in the instruction 
which was replaced at Step 08 in Step 16, and the three missing instructions 
as 17, 18, and 19. Then key in GTO 09 as Step 20. 
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CALCULATOR PROBLEMS 

!• Using black and white Infrared film and a suitable filter one can 

photograph within the wavelength range 500-650 nm. If you were to attempt 

to photograph the reindeer in problem 3 (main text), how much power would 

2 

be reflected in these wavelengths by a m of: 

a) The reindeer horn at T = 32**C? 

b) The outside of the reindeer's fur at T = 5'*C? 

c) The external environment at T = 5**C? 

At present. Infrared film is not "fast" enough to record thermal 
luminescence of an object at temperatures below 250**C, 

d) How much power would be available in this wavelength range 

from a black body at this temperature? 

The solution to this problem will require the numerical integration of 

f6.5xiO"^ 



I = 3.7A X 10"^^ 



-5[e ^^/^^T . l]-^dX 



5x10"^ 



which can be derived in a manner analogous to the derivation of Eq. 1 but 
beginning with the solution to problem 1 (main text), that is 

2. The sensitivity of the compound eye of a honeybee (worker) extends 
further into the ultraviolet than does the sensitivity of the human eye. 
The bee eye is sensitive from about 300-650 nm, and the human eye is sensitive 
from about 400-700 nm. \Jhat is the difference in the power available from the 
sun in the wavelengths perceptible to these two organisms? 
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The solution to this problem will require the numerical .iitegration of: 



I = 3.74 X lO"-"-^ 



X -5 ^^2.52x10-^2^^-1 



with the limits of integration chosen to c'^incide with the spectral ranges 
provided above. 

• • 

Mhat is the biological significance of the greater ultraviolet sensitivity 
of the bee eye? 

3. Photosynthetic bacteria differ from most algae and higher plants 
in thcit they possess an auxiliary "electron pump" in their photosynthetic 
' . ' which allows them to utilize radiation of longer wavelengths for 
p.' ^-ntnesis. The absorption spectra for the photosynthetic organs of 
these organisms reveal the following spectral ranges: 

a) 400-700 nm green plants 

b) 780-840 nm ' 

c) 840-910 nm - various purple bacteria 

d) 970-1050 nm^ 

Using the integral in problem 2 above, calculate the power available 
from the sun in each of these wavelength ranges. 

For an enlightening discussion of the evolutionary significance of the 
wavelength ranges utilized by these bacteria, see Bjorn 1976 (p. 235). 
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CALCULATOR PROBLEM SOLUTIONS 



!• The Integral to be evaluated Is 

.-7 



I « 3.74 X lO'*"'-^ 



r6.5xi0 
-7 



5X10 

For the four parts of this problem the value of T differs, which means 
that the exponent of the exponential of the integral will be different for 
each part. 

a) T = 305^K; ^ = 4.718 x lO"^ 

b) T = 278^K; = 5.1757 x lO"^ 

c) T « 268^K; 3^3^33 ^ ^q-5 

d) T = 523^K; 2.751 x 10*"^. 
For the SR-52: 

As explained in section C of Appendix 2. II, one must alter the program 
in Table A-3 at step 095 which is the function defining subroutine. A suitable 
program is given in Table A-5. The preloaded registers in Table A-5 are 01, 
02, and 07. Register 07 requires the values of hc/kT given above. Note that 
the va;:.je of the integral calculated by the program must be multiplied by 
3.74 X 10 to give the correct value for the power. 
For the HP-25: 

As explained in section C of this appendix, one must alter the program in 
Table A-4 at step 12. A suitable program is given in Table A-6. The preloaded 
registers are 1, 2 and 7. Register 7 requires the values of hc/kT given above. 



Note that the value of the integral calculated by the program must be multiplied 

8i 



by 3.74 X 10 to give the correct vaiue.f or the power. 
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The values for the Integral that I obtained were 



a) 


1 = 


2.A1 


X 


10-9 


b) 


I = 


2.56 


X 


10-^2 


c) 


1 = 


9.A5 


X 


10-^^ 


d) 


I = 


6.01 


X 


10-^ 
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Table (SR-52) program Changes to Table A'3 



00 



01 



02 



ERIC 



.5.7X10"^ 6.5x10"''/ X 



03 

i 



04 

Ax 





(>reloaded 






Cocl^ 


Key 


095 




sio 


096 




1 


097 


00 


9 


098 


4^ 




099 




5 


lOO 


6? 


X 


101 


53 


( 


102 


53 


( 


103 


43 




104 


00 


0 


105 


0/ 


F 


106 . 


55 


f 


107 


43 


RC'L 


108 


01 


1 


109 


09 


9 


110 


54 


) 


111 


22 




112 


23 




113 


75 




114 


01 


1 


115 


54 


) 



05 



06 



i-1 



07 

kT 
+ 



08-18 



19 

used for 
f(x) 



Comments 



Store 



^1 



Register 19 



calculate \^ 



Rec?ll hc/kT 



R<^.calJ X.' 

1 



Calculate hc/XkT 



Calculate exp[hc/XkT] - 1 



Calculate X [exp hc/XkT>.i] 
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00 01 02 03 04 05 06 07 08-18 19 

^ ^ ^i-1 w — ^""^ ^^"^ 

preloaded preloaded 



S tep Code Key Entry Comments Comments 

116 95 « Calculate X'^texp hc/XkT-1] - 1 

117 20 1/x Calculate X_j^"^[exp(hc/kT)-l]"'"'' 



118 56 rtn 



o 

ERIC 



I 

I 









Table 


A-6 (IIP-25) 


Program Chaugt-f to Table A-4 












Register Contents 




0 




1 


2 


3 


.4 


5 6 


7 






a b 

preloaded 


^i 




i-1 


ill-/ l\iX 


Step. 
12 


Code 


Kev Entry X 


Y 


Z T 


Comments 




31 


ENTER 




h . 






13 




31 


ENTER 


\ 


h 






14 




05 


5 


5 


h 


\ 




" ■: 


14 


03 


fy 






h 


Calculate X. 


It 
17 




21 
07 


s < y 
RCL 7 


hc/kT 


, 5 
^i 

^i 


\ 5 


Recall hc/kT 






7.1 


X < y 


i 


hc/kT 




Calculate hc/XkT 


19 
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hc/XkT 


X. 


,5 
\ 


Calculate hc/XkT 


20 


15 


07 


X 


exp[hc/XkT3 






Calciilate expihc/Xk 


21 




01 


1 


1 


exD F 1 


X.^ 




22 




41 




exp[]-l 


X 




Calculate expthc/Ak 


23 
24 


15 


61 
22 


X 

1 


[exp[~]]"^ 






he 

Calculate [exp Cj^]^) 


25 


24 


04 


RCL 4 


Ax 


f (Xj^ + -y) 




in Table A-4 






Continue k ■ -li 


steps 26-33 


as in steps 


23-30 in Table A-4 








For 


this program, Step 31 Table A-4 must be changed to: 






13 


38 GTO 38. 
Continue T-rith 


steps 35-49 


as in steps 


32-46 in Table A-4. 
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2c The integral to be evaluated is the same as in problem 1, Thus one 

can use the programs given in Table A-5 (or A-6) preloading register 07 
(or 7) with: 

« 2.52 X 10-^^ 
The value for the integral that I obtained was 

a. 300-650 nm: I = 6.38 x 10^^ 

b. AOO-700 nm: I 5.85 x 10^^ 

3. The integral to be evaliiated is the same as in problem 2. The values 
that I obtained for the different ranges of integration were 

a. 400-700 nra; I = 5.85 x lO^^ 

b. 780- 810 nm: I = 8.02 x lO^""" 

c. 840- 910 nm: I = 8.123 x 10^^ 

d. 970-1050 nm: I = 6.85 x lO^-"" 
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Appendix 3: SOLUTIONS TO PROBLEMS 

Problem 1 

Planck's formula in terms of \: 
E(v) = 8Trhv^/c^ (e hv/kT - 1)"^ . 
So ECA)dX = -E(v)dv Oninus sign because dX and dv are of 

opposite sign) . 

Since V = c/X for all wave phenomena. 



dv/dX = -c/X^. 



E(X) = ECv) dv/dX = E(v) e/X^ 



= c/x2 



87rhv 



, 3 
^ c 



4i 



[e hv/kT - I]"'"'- 
= (87ThcA^)[e^^/XT - l]*""^. 

Problem 2 

The general formula is 

X T = b or X = b/i\ So 
m m 



^ 1 2,8978 x 10""\ k c no •.n-7 
a) X = = 5.08 X 10 m 



5.7 X 10^ K 



b) Since 37°C = 300 K 



, _ 2.8978 X 10"\ K _ q .cq ^ m"^ 
A = =; = 9.659 X 10 m 

° 3 X 10^ K 



c) Since lO'C = 283 K 

_ 2.8973 X 10~\ K _ , ^9/ v in-3 
X = ' 5 = 1.024 X 10 m 

^ 2.83 X 10 K 



90 



ERIC 



82 



d) Since -33''C = 240 K, 



^ 2.898 X io"\ K ^ , ^ ,^-3 
= 7i = 1.207 X 10 m. 

2.4 X 10^ K 



Problem 3 



As In example 2 

a) S = 0|t^^ - T^^J = 5.67 X I0~^j^(273 + 32)^ - (273 - 10)'^ W m" 
= 219.4 w m"^ 



H = E(SA) = 219.4 W m~^(.l m^) = 21.9 



b) E = a^lj" - T^^J = 5.67 X 10"^jj:305)^ - (228)^ w m~^ 
= 337.4 Win-2 

H = E(SA) = 337.4(7.5 x 10~^) = 25.3 W. 

Problem 4 . 

4 8 

= aT^ ^/v 2 = 5.67 x 10"^^5.70 x 10^) iZJ^JLlOf. 
^ ^ ^ (1.49x10^^)2 

= 1.359 X lO'' W m''^. 

Converting to calories 

Eg « 1.359 X 10^ W ia"2=: 1.359 Wm''^ |^.0239^^^/j[ lO^V/ 
-1 -2 

= 1.95 cal min cm . 



Problem 5 



As In Example 4: 
a = 0.61 X/r 
A = 5.08 X lO^'^m 

-4 

r 1.0 mm = 1.0 x 10 m 
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-4 

So a = 3,17 X 10 rad. 



Again as in Ex. 4, 

6x = ad « (3.17 X lo"*^ rad) (17 ram) 
-3 

6x « 5.4 X 10 mm 
= 5.4 X lo""^ m. 



Problem 6 

From the geometry of the situation it is obvious that 
d = aR (i) 

^e a is the angle of the cone which is the onomatadia. If this angle 
is at the minimum given by Rayleigh's criterion: 

a = 0.61 X/r = 1.22 X/d (ii) 
Combining equations Ci) and Cii) and solving for d, one gets 

d = (1.22RX)^/2 . 
Substituting in R = 1 mm and X = 400 nm gives 
d = 2.2 X 10"\. 

Problem 7 

X = 9.67 X 10"^m 

r = striking distance ^ 1 m 

s = accuracy of strike Cr 2 cm 
-2 

(J) = s/r « 2 X 10 rad. 
Using the Rayleigh criterion, the separation of the pits, d, is given by 

d = 2.(0.61)X/a = 1.22 ^'^^ 1° « 6 x 10^^^ 

2 X lO"" rad 

^ 6 cm. 

This theoretical prediction corresponds very closely to the actual 
observed separation of the infrared sensing organs of a pit viper. 
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Problem 8 



The echo will return after: 

T = d/c^ = lOm/350 u s""*" - 28.6 ins. . 
The insect detects a frequency: 

which is reflected to the bat and detected at the frequency: 
^L-«3.V,,V(C3-V,,.mi|j:3|0i|] 

= 51.74 kHz. 
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